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Abstract 

b : 

We reconsider the 187 Re bound on the variation of the fine-structure constant. We 
combine the meteoritic measurement with several present-day lab measurements to firmly 
establish the quantitative upper limit to the time variation over the age of the solar system. 
We find that the relative change of the fine-structure constant between its present value and 
a of ~ 4.6 Gyr ago is consistent with zero, Aa/a = [—8 ± 8 (ler)] x 10 -7 . We also rederive 
this bound in models where all gauge and Yukawa couplings vary in an interdependent 
manner, as would be expected in unified theories. 



The prospect that the fundamental constants of nature vary in time is exciting and rep- 
resents a radical departure from standard model physics. Recent attention to this problem 
has focused on the indication that the fine structure constant was smaller at cosmological 
redshifts z = 0.5-3.5 as suggested by observations of quasar absorption systems [l]-[4]. 
The positive result of Aa/a = (0.54 ± 0.12) x 10~ 5 , where Aa is defined as the present 
value minus the past one, is based on the many-multiplet method which makes use of the 
a dependence of the relativistic corrections to atomic energy levels. 

However, there exist various sensitive experimental checks that constrain the variation 
of coupling constants (see e.g., [5, 6]). Among the most stringent of these is the bound 
on | Aa/a | extracted from the analysis of isotopic abundances associated with the Oklo 
phenomenon [7]-[10], a natural nuclear fission reactor that occurred about 1.8 billion years 
ago. While the Oklo bound |Aa/o;| < 10~ 7 is considerably tighter than the "observed" 
variation, the Oklo phenomenon occurred at a time period corresponding to z ~ 0.14. It 
is quite possible that while a varied at higher redshifts, it has not varied recently. That is, 
there is no reason for the variation to be constant in time. Big bang nucleosynthesis also 
provides limits on Aa/a [11, 12, 13]. Although these limits are weaker, they are valid over 
significantly longer timescales. 

Bounds on the variation of the fundamental couplings can also be obtained from our 
knowledge of the lifetimes of certain long-lived nuclei. In particular, it is possible to use 
relatively precise meteoritic data to constrain nuclear decay rates back to the time of solar 
system formation (about 4.6 Gyr ago). Thus, we can derive a constraint on possible vari- 
ations at a redshift z ~ 0.45 bordering the range (z = 0.5-3.5) over which such variations 
are claimed to be observed. The pioneering study on the effect of variations of fundamental 
constants on radioactive decay lifetimes was performed by Peebles and Dicke [14] and by 
Dyson [15]. The isotopes which are most sensitive to changes in the Q value are typi- 
cally those with the lowest /5-decay Q-value, Qp. The isotope with the smallest Qp value 
(2.66 ± 0.02 keV) is 187 Re. 

In Dyson's analysis [15], determination of 187 Re lifetime was based on (1) isotopic mea- 
surements of molybdenite ores [16], (2) isotopic measurements of iron meteorites [17], and 
(3) a direct measurement of the decay rate [18]. The decay constant was found to be 
A = (1.6 ±0.2) x 10" 11 , (1.4 ±0.3) x 10~ n , and (1.1 ±0.1) x lO" 11 yr^ 1 , respectively. The 
discrepancy between the direct measurement and the geophysical ones led to speculations 
that this could be due to a time variation in the fine structure constant. Note that this 
would indicate a decreasing value of a from the past to the present, i.e., 5a < 0. However, 
the direct measurement [18] was complicated by many technical issues, and a later labo- 
ratory measurement [19] gave A = (1.5 ± 0.2) x 10 -11 yr" 1 , which was consistent with the 
geophysical measurements. As a conservative upper limit, Dyson [15] concluded from the 
above measurements that \Aa/a\ < 5 x 10~ 6 or \a/a\ < 5 x 10~ 15 yr -1 over the last Gyr. 
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This was less stringent than the Oklo bound [8]. 

Recently, we have reconsidered the constraint on time variations of the fundamental 
couplings based on meteoritic measurements of the lifetimes of radioactive isotopes [20] . In 
particular, we made use of the dramatic improvement in meteoritic analyses of the 187 Re- 
187 0s system in iron meteorites. The present abundances of 187 Re and 187 0s in a meteorite 
are 

(i 87 Re) = ( 187 Re) ie xp(-At), (1) 

( 187 0s) = ( 187 0s)j+( 187 Re)j[1 _ exp( _ At)]) (2) 

where the subscript "i" denotes the initial abundance and t is the age of the meteorite. 
Equations (1) and (2) give 

(i87 0s ) = ( 187 0s), + ( 187 Re)[exp(At) - 1]. (3) 

Thus, there is a linear correlation (called an isochron) between the present abundances of 
187 0s and 187 Re (measured relative to the reference stable isotope 188 0s, which does not 
receive any decay contributions). 

The slope of this correlation, exp(At) — 1, can yield the 187 Re decay constant A if the 
age t is independently known. Radiometric data support the idea that the Group IIIA iron 
meteorites were formed at about the same time (±5 Myr) as the angrite meteorites [21], 
which have a precisely determined 207 Pb - 206 Pb age of 4.558 Gyr [22]. The combination of 
this age with the slope of the well defined 187 Re - 187 0s isochron for IIIA iron meteorites gave 
A = (1.666 ± 0.009) x 10~ n yr" 1 [21] (see also [23, 24]). The dominant contribution to the 
la uncertainty is associated with the non-stochiometry of the Os salt used for calibration 
(~ 0.5%). Minor uncertainties include the Re - Os isochron slope, the 235 U and 238 U decay 
constants [25], the time difference between closure ages for U - Pb pair in angrites and the 
Re - Os pair in IIIAB iron meteorites. 

The /5-decay of 187 Re is a unique first fordidden transition, for which the energy depen- 
dence of the decay rate can be approximated as [26] 

A oc G\Qlm\. (4) 

As we showed previously [20], considering only the variation of the Coulomb term in Qp, 
we have 

AA =3 Aa £ ^3(20MeVWAa^ 2xl04 /A<»X 
A Qp \ Qp J V a ) \ a J 

Thus, given a determination of the possible variation in A, we can obtain a limit on the 
variation of a. In our previous work [20], we assumed that the variation of A cannot exceed 
the accuracy of the meteoritic measurements, AA/A < 0.5%. This gave Aa/a < 3 x 10~ 7 
or a/a < 6 x 10~ 17 yr _1 over a period of 4.6 Gyr, assuming a linear evolution with time. 
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Two issues concerning the above limit requires attention. The meteoritic ages used to 
derive the decay constant of 187 Re are determined in part by the decay constant of another 
radioactive isotope 238 U, which also depends on a. However, as noted earlier [20], this does 
not change our limit on AA/A because the decay constant of 187 Re is far more (by a factor of 
40) sensitive to changes in a than that of 238 U. On the other hand, the assumption that the 
variation of A cannot exceed the accuracy of the meteoritic measurements is not justified. 
Strictly speaking, this assumption is valid if measurements of comparable precision are 
made for meteorites with a large spread (e.g., several Gyr) in age and yield consistent 
decay constant for 187 Re. In practice, accurate measurements are available only for iron 
meteorites with essentially the same age. Thus, our earlier limit based on AA/A < 0.5% 
over a period of 4.6 Gyr was over-restrictive. 

To obtain a solid limit on AA/A, we can compare the meteoritic measurements of A, 
which cover 4.6 Gyr, with direct laboratory measurements, which give the present value 
A(to) = Ao- Indeed, there is a reasonably accurate direct measurement which yields Ao = 
(1.639 ± 0.025) x 10~ n yr" 1 [27] (see also [28]). The variation in A is then AA = A - A = 
(-0.027 ±0.026) x lO" 11 yr" 1 , or 

AA 

— = -0.016 ±0.016. (6) 
A 

As one can see, there is an insignificant trend in the negative direction. However, it is worth 
mentioning that the implication of the quasar absorption system measurements would be a 
positive AA, under the natural assumption that a(t) is a smooth monotonic function. This 
result and Eq. (5) give Aa/a = [—8 ± 8 (ler)] x 10~ 7 , which yields a 2a limit 

-24 x 10~ 7 < — < 8 x 10~ 7 . (7) 

a 

This is weaker than the limit reported in [20] by a factor of ~ 3 for a positive AA and by 
a factor of ~ 8 for negative AA. Even so, the O(10~ 6 ) limits imposed by the meteoritic 
data still provide strong constraints on models of a(t) and severely restrict possibilities to 
accomodate the claimed variation of a based on observations of quasar absorption systems 
(see, e.g. Ref. [29]). This result is also about 6 times stronger than that given in [28] 
based on the meteoritic data of [24]. (We note that a more recent measurement [26] gave 
Ao = 1.68 x 10~ n yr -1 with a larger systematic uncertainty of ~ 3%. Using a weighted 
average of the two direct measurements of Ao does not affect the above result significantly.) 

It is important to note that this approach, i.e. the comparison of the meteoritic mea- 
surement of A with Ao, will benefit from any improvement in the determination of Ao which 
is required for the application of the Re - Os geochronology [30]. On the other hand, it is 
also important to recognize the limitations of this limit. Strictly speaking, if A varies in 
time, then the meteoritic measurement really represents the time average of X(t) over the 
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age of the meteorite. As such, the redshift at which the limit can be applied will in principle 
depend on the specific functional dependence of \(t). For example, if we consider a simple 
mononomial form, \(t) — Ao ~ {to — t) n , then the limits (6) and (7) should be applied at a 
reduced look-back time of (to — t) = 4.6Gyr/(n + I) 1 /". As noted above, a look-back time 
of 4.6 Gyr corresponds to a redshift of z ~ 0.45 for the assumed cosmological model, hence 
the limit would be applied (in a model dependent way) at a somewhat smaller redshift. 
For example, in the linear case, n — 1, one can either apply our limits above at a time 
t — t — 2.3 Gyr, or one can apply a relaxed limit (by a factor of 2) at a look-back time of 
4.6 Gyr. Because of the model dependence, it may be possible with a careful choice of a 
time dependence for A to obviate this limit [31] due to the limit being tied to a time average 
of the decay rate. 

As we discussed previously, in the context of unified or string-inspired theories, all 
Yukawa couplings and gauge couplings depend on the same moduli fields and the change 
in the fine structure constant typically implies a change in other couplings and mass scales 
[12, 32]. The dominant effects are found in induced variations of the QCD scale A and 
the Higgs expectation value v. The variations ^ ~ 30^- and — ~ 80^ [12, 33, 34] are 
translated into variations in all low energy particle masses. In short, once we allow a to 
vary, virtually all masses and couplings are expected to vary as well, typically much more 
strongly than the variation induced by the Coulomb interaction alone. In Ref. [20], we 
adopted |AA/A — Av/v\ ~ 50Aa/a;. The contributions to Qp from (the kinetic energy) T 
and (the nuclear potential energy) V, which scale with A, are comparable to that from the 
dominant Coulomb term (C), which scales as aA. As changes in A are O(30) times larger 
than that in a, we can estimate 

A Q p v Q/3 A A 



which gives 



-8 x 10~ 8 < — < 3 x 10" 8 , (9) 
a 



over a period of 4.6 Gyr. 

In summary, we have revisited the bounds on the change of the fine-structure constant 
from the 187 Re meteoritic measurements. We have shown that these limits are determined 
by the accuracy of the laboratory determination of the decay rate. This results in a bound, 
Aa/a = (—8 ± 8) x 10~ 7 going back to ~ 4.6 Gyr ago, that is of significant importance for 
all models of a(t). 
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